Dynamical nuclear spin polarization and the Zamboni effect in gated double quantum dots 
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A dynamical nuclear polarization scheme is studied in gated double dots. We demonstrate that a small polar- 
ization (~ 0.5%) is sufficient to enhance the singlet decay time by two orders of magnitude. This enhancement 
is attributed to an equilibration process between the nuclear reservoirs in the two dots accompanied by reduced 
fluctuations in the Overhauser fields, that are mediated by the electron-nuclear spin hyperfine interaction. 
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Electron spins localized in semiconductor quantum dots 
have been intensively investigated in recent years due to 
their potential use in quantum information processing and 
spintronics. 1 Several experimental 2 ' 3,4 and theoretical 5 ' 6 ' 7 ' 8 ' 9,10 
studies have identified the hyperfine (HF) interaction between 
an electron spin and the surrounding nuclear spins as one of 
the main sources for electron spin decoherence in low tem- 
perature GaAs quantum dots, leading to T 2 * on the order of 
10 — 25ns and T2 on the order of [is. These values are or- 
ders of magnitude shorter than the spin relaxation time, which 
approaches tens of milliseconds in these systems. 2 

Several strategies have been suggested to alleviate elec- 
tron spin decoherence via nuclear spins, including spin- 
echo techniques (to remove inhomogeneous broadenin g 4 ' 7 - 10 ), 
nuclear spin state measurement (to narrow the Overhauser 
field distribution, 1 ') and nuclear polarization (to reduce phase 
space for nuclear spin dynamics. 12 ) Nuclear spin polarization 
is also valuable for state initialization in NMR quantum com- 
puting and for utilizing collective nuclear states as long-lived 
quantum memory. 13 So far, optical pumping has produced up 
to ~ 60% nuclear polarization 14 in interface fluctuation GaAs 
dots, while spin transfer via hyperfine mediated spin-flip scat- 
tering in the spin-blockade regime in gated GaAs dots has 
led to ~ 1% polarization. 15 ' 16 Among the limiting factors in 
these dynamical polarization schemes are the large difference 
in the Zeeman energies of the electrons and nuclei making 
joint spin flip processes energetically unfavorable, and nuclear 
spin diffusion due to dipolar interaction. In contrast, theoreti- 
cal studies have shown that in order to achieve a sizeable en- 
hancement of single electron spin decoherence time via phase- 
space squeezing, a nuclear polarization of more than 99% is 
required . 8 ' 9 ' 17 

In this paper we show that HF interaction can be exploited 
to dynamically polarize the nuclear spins in gated double dots. 
Most interestingly the relaxation time of the two-electron spin 
singlet state dramatically increases without the need for nearly 
complete nuclear polarization. This suppression of relaxation 
is achieved by an equilibration process in the nuclear reser- 
voirs in the two dots and a reduction in the fluctuations of 
their Overhauser fields. We have dubbed this effect, mediated 
by the HF interaction during the polarization cycles, as the 
nuclear Zamboni effect^ 

We study the dynamics of a system of two electrons local- 
ized in a gated double dot interacting with two nuclear spin 
baths within the framework of the Hamiltonian 



where we have neglected nuclear-nuclear dipolar coupling in 
the current study. For the orbital part we adapt the Hund- 
Mulliken approach 12 ' 19 to solve for the electronic states in 
the gated dot configuration, where iT rb includes the single- 
particle Hamiltonian and the Coulomb interaction. The rel- 
evant Hilbert space of 7J or b is spanned by four two-particle 
states, {5(2,0), S(Q, 2), 5(1, 1), T(l, 1)}, which consist of 
the separated singlet and triplet and the two doubly occu- 
pied singlet states. Indices indicate the number of con- 
fined electrons in the (left,right) dot. We neglect the dou- 
bly occupied triplet states as their energy is much higher for 
the structures under study. 2 - 4 The Zeeman interaction Hz — 
• Yli=L R&i> Wltn 9 = — 0.44 and \ib being the 
Bohr magneton, splits the triplet states. The resulting energy 
diagram and exchange interaction J near the (1,1) to (0,2) 
charge transition are shown in Figs. [Tk,b for B = 100 mT 
(Ez = 2.5/ieV), half interdot distance a — 1.9as, and dot 
confinement ojo = 120/icV, corresponding to the experimen- 
tal parameters in ref.0. 

For the HF interaction we consider the contact term be- 
tween the two electrons and the surrounding nuclei 
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where is the HF coupling constant with the fcth nucleus 
in the ith dot. Here h = | (h L + h R ), Sh = ~ (h L - h R ) 
and S = Sl + Sr, 8S = Sl — S R are the sums and dif- 
ferences of the nuclear fields and the electron spins in the two 
dots, respectively, where hj = J^k ^i ^i ls me nuclear field in 
dot i. We assume l£ — ^ for simplicity. The Hamiltonian in 
Eq. ([T]i conserves the total spin and can be block diagonalized 
in each of the Hilbert subspaces defined by the eigenvalues of 
the operator J z = 5£ + S R + J2i k ^ ■ ^ n or der to make 
the numerical effort for a larger number of spins tractable, we 
assume a constant HF coupling for all the nuclei in each dot, 
allowing us to write the Overhauser fields as = (7j/iV)Ij, 
i = L,R, where 7i = J2k A i ~ 100 ^eV is the total HF 
coupling, N is the number of nuclei per dot, and Ij is the col- 
lective spin operator for dot i. This approximation provides 
two more integrals of motion, namely the two SU(2) Casimir 
operators I|, I R , and enables us to further divide the Hilbert 
space, making the complexity of the problem scale polynomi- 
ally with N instead of exponentially. 

We have tested the validity of the uniform HF coupling 
approximation for the nuclear polarization dynamics stud- 
ied here by slicing the dot into concentric rings. Assigning 
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different HF coupling for each ring, and assuming no inter- 
ring dynamics, we find that the dynamics described below 
are largely unaffected by this averaging procedure. 20 Studies 
of nuclear spin diffusion in a quantum dot also verify that 
inter-ring nuclear spin dynamics is significantly suppressed 
for smaller quantum dots. 21 Our approach enables us to study 
the interplay between HF and exchange effects within a uni- 
fied theory through exact diagonalization of the Hamiltonian 
([T]i. The system dynamics under any applied gate-pulse are 
calculated without resorting to the quasi-static approximation 
that may not be appropriate to describe nuclear polarization 
dynamics^ 2 - At the same time, our use of collective spin states 
enables us to consider a substantially larger number of nu- 
clei (~ 1000 per dot) as compared with previous numerical 
studies, 6,7 which is important in determining the scaling prop- 
erties of the dynamical features with N. 

The nuclear spin configuration in each dot is represented 
in the basis of collective Dicke states \j, m), where < j < 
N/2 is the total spin of the state (or cooperative number) and 
\m\ < j is the z projection of the total spin. The initial nuclear 
spin configuration in the two dots is written as 

N/2 N/2 
|V(0)>nuc= £ W(I L )\I L ,II°)®J2 W(I R )\I R ,F R °), (3) 
Il=0 Ir=0 

and will be denoted henceforth as (T£°, I R ) W where If in- 
dicates the initial polarization of the nuclear configuration in 
dot i. The distribution weights w(Ii) are assigned statistically 
by the number of possible ways in which one can arrange N 
spins into singlets, triplets, quintets, etc. We find that these 
weights obey a Gaussian distribution peaked at T = \/N/2 
whose width is v2N, as shown in Fig.QJi 2 ^ 

An example of the proposed dynamical polarization pulse 
sequence is shown in the lower panel of Fig. QJ. 16 The key is 
to drive the system through the 5 — T + resonance with differ- 
ent speed in the two directions, so that only in one direction 
can the electron-nuclear spin flip occur. An 5(0, 2) state is 
prepared by positively detuning the double dot (P) to enabe 
electron exchange with the leads, whose Fermi level is above 
the 5(0, 2) but below the doubly occupied triplets. The elec- 
trons are then separated using rapid adiabatic passage, where 
the bias is swept to a negative detuning quickly relative to the 
HF coupling but slowly as compared to the electron tunneling 
between the hybridized singlet states. Our simulations show 
that the adiabaticity requirement is met with sweep times of 5 
ns. The bias is then swept back slowly through the 5 — T + an- 
ticrossing (D, which can be identified experimentally 4 ) where 
the HF interaction mediates electron spin flip flop with the nu- 
clear spin baths. Choosing point 5' to be far enough from the 
5 — To degeneracy, the electron spin state is always flipped 
from a singlet to a T + , thus polarizing the nuclear spin baths. 
Finally, the system is swept back to point P where the triplet 
state relaxes quickly through electron exchange with the leads, 
and a new 5(0, 2) state is prepared for the next cycle. The 
singlet preparation at the end of each cycle is simulated by 
partially tracing the electronic subsystem and applying a di- 
rect product of the resulting nuclear density matrix with the 
electronic configuration at point P. To study effects of the po- 




Bias parameter 

FIG. 1: (Color online) (a) Upper panel: Orbital energy diagram for 
the double dot near the (l,l)-(0,2) transition vs. a bias parameter, 
proportional to the inter-dot bias gate potential. Shown are the hy- 
bridized singlet states (black curves), and split (1,1) triplet states T_ 
(dash-dotted green), To (blue), and T+ (dashed red). • denotes the 
S — T+ degeneracy point. Lower panel: Polarize, separate and mea- 
sure pulse sequence. Three polarization cycles are shown followed 
by a separation time and a measurement of the singlet probability. 
The letters at the bottom indicate biasing points discussed in the main 
text, (b) Exchange energy as a function of bias, (c) Center location 
(blue) and FWHM (red) of the nuclea r spin state distribution as a 
function of N. The dashed lines are \J N/2 and V2N, correspond- 
ing to the center location and FWHM of the distribution, respectively. 



larization cycles on the relaxation of the electron singlet state, 
we add a measurement cycle where the dots are negatively 
detuned to the 5 — To degeneracy (5) for a separation time 
rg, followed by a measurement of the singlet probability Pg 
(t~m ~ 5/xs)!^ 4 - To enhance the efficiency of the polarization 
process, we perform a nonlinear bias sweep, spending a sub- 
stantial part of the cycle in the vicinity of D (see Fig.QJ). 

Before presenting the simulation results for the polariza- 
tion scheme, we discuss the envisaged impact of this pro- 
cedure on the decay time of the electron singlet correla- 
tions. Figure [2] depicts the time evolution of Pg when 
the electrons are prepared in a singlet state and placed at 
point 5, for several values of N. In Fig. [2^ the initial 
nuclear state is (Q, — mutjViV, Ir]) w , representing a sta- 
tistical polarization difference between the two dots. The 
time axes are multiplied by sj TV/10 5 , indicating that the 
decay time scales like 1/y/N. The decay time of 25 ns 
agrees well with the experimental findings, 4 ' 24 and its scal- 
ing with N corresponds to the decoherence time behavior 
found in ref. 8, since we scale the HF coupling constant with 
a fixed number of nuclei (10 5 ) rather then with N. Sim- 
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ilar dynamics is found for other initial nuclear states, e.g., 
(— min[viV, Ii], — min[2v A /V, Ir])w Fig. Eh shows the time 
evolution of Ps for the fully polarized nuclear configuration 
(—II, —Ir)w In this case Ps decay times are two orders of 
magnitude longer reaching 3/is and they do not scale with TV. 
Other equally polarized configurations do not always present 
an enhanced singlet coherence, indicating that the fully polar- 
ized nuclear state is characterized by a narrower distribution in 
addition to having 8h z = 0. We stress that the term fully po- 
larized does not suggest that all the spins are polarized, since 
one is limited by the total spin of each collective state within 
the weighted distribution. In fact, the total attainable polar- 
ization is p = 1.4533/ vN reaching a value of 0.46% for 
TV = 10 5 , which is consistent with recent experiments.-^^ 
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FIG. 2: (Color online) Electron singlet probability as a function of 
separation time for TV values ranging from 100 to 1296, and J = 0. 
(a) Initial nuclear configuration (0, — min[yTV, Ir])w Time axes 
are scaled like y/ N/10 5 , so that the presented decay times corre- 
spond to TV = 10 °. The rise in Ps indicates decayed oscillations, 
which are an artifact of the uniform coupling approximation and can 
be eliminated by performing dynamics averaging over different HF 
couplings, (b) Initial nuclear configuration (—II, —Ir)w 

The short time dynamics can be understood within a 2 x 2 
effective Hamiltonian for the S z — subspace, H c g — 
4(1 + t z ) + 5h z r x , where r is the pseudospin operator 
(\S) -► \t z = -1), |T > - \t z = l)). 8 Even when (Sh z ) = 
as for the fully polarized state in Fig. [2b, the spin dynamics 
do not vanish altogether. The reason is that applying H c h 
on each of the collective states in the weighted distribution re- 
sults in different eigenvalues and their dynamics do not cancel 
out, giving rise to quantum fluctuations. In both cases shown 
in Fig. [2] there is a long time (~ 2.5/xs for B = 100 mT) 
envelope decay attributed to higher-order corrections to H c g 
that are contained in Eq. ((TJ. This envelope scales like 1/TV 
and could therefore govern the dynamics of the fully polar- 
ized state in the large TV limit. Its long timescale is a con- 
sequence of the large Zeeman splitting as compared with the 
Overhauser fields and it can be made longer using a larger 
magnetic field. For the TV values we are considering, these 
corrections are only observed for the state (0, 0)^ for which 
all other dynamics are shut down. 

Our results also agree well with the experimental findings 24 
and analytical results 8 for J > 0. These include preservation 



of the singlet correlations over a long timescale in the limit of 
J ^> E nuc , and the appearance of damped oscillations in Ps 
with a saturation value that depends on E nuc / J in the interme- 
diate regime J ~ E nuc . We also find the long time (rg ^> T 2 *) 
value of P s to be 1/2 for B > B nuc and 1/3 for B < B nuc , 
in agreement with semiclassical theory 5 . 

Now we investigate the effects of the polarization cycles on 
the electron spin states by separating the electrons every four 
cycles to calculate Ps as the polarization progresses. Bias 
changes require calculating the evolution separately for each 
bias, using the resulting density matrix at each step as the ini- 
tial condition for the subsequent step. The numerical effort is 
thus much more demanding and we are limited to several tens 
of nuclei per dot. Figs.[3h-c show Ps calculated for the initial 
nuclear configuration, after 20 polarization cycles, and after 
100 cycles. The singlet decay times show a gradual enhance- 
ment as the nuclear polarization builds up, and their scaling 
with TV gradually shifts from 1/-//V (Fig. Hi) to 1 (Fig. Hfc). 
An enhancement of factor 300 is obtained for the singlet de- 
cay times that reach ~ 8/is when the polarization process is 
complete. The corresponding nuclear polarizations in the two 
dots, shown in Fig.fJJl, equilibrate during the polarization pro- 
cess. This equilibration effect is robust to any of our choices 
of initial nuclear configuration, and the degree of equilibration 
depends on the symmetry of the double dot. The equilibration 
between the two nuclear spin configurations and the narrow 
distribution of the Overhauser fields formed during the polar- 
ization process are responsible for the prolonged singlet decay 
time. 
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FIG. 3: (Color online) Singlet probability vs. separation time, calcu- 
lated after performing polarization cycles, for several TV values. The 
initial nuclear configuration is (0, — minfvTVi Ir])w and J = 0. (a) 
Initial Ps without polarization. Time axes are scaled like ^/TV/10 5 
(b) Ps calculated after 20 polarization cycles. Time axes are scaled 
like (TV/10 5 ) ' 28 (c) P s calculated after 100 polarization cycles. 
Time axes do not scale with TV (d) The corresponding nuclear po- 
larizations in the left (L-red lines) and right (R-blue lines) dot. The 
vertical dotted line corresponds to the time elapsed after 20 cycles, 
at which Ps shown in (b) was calculated. 



Dipolar interaction between the nuclear spins can break 
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the weighted distribution of the spin collective states, and 
compete with the equilibration process. The number of po- 
larization cycles needed to complete the process is n cyc = 

^Zifio ~ 1.4533V0V. The spin transfer time at point 
D is found to be 15.5/Js/y/N. As long as the S — T + degen- 
eracy can be determined accurately, much of the cycle time 
in the large N limit is spent to prepare a singlet at P. Tak- 
ing tp = 200 ns for the electron exchange with the leads, the 
resulting time to complete the polarization is ~ 120/is. 

A simple phase space argument seems to suggest that the 
interdot Zamboni effect should occur naturally through nu- 
clear spin exchange between the two baths mediated, e.g., by 
one or two electrons in the double dot. While the fully polar- 
ized state has a single configuration, the non-polarized state 
is highly degenerate having C^ 2 configurations (assuming 
equal HF coupling for all the nuclei). Even for states with 
statistically similar dot polarizations the phase space differ- 
ence could drive an equilibrating process. For example, the 
ratio between the number of configurations in (0, VN) W and 
(VW/2, VN/2) w is e" 1 / 4 which, as our simulations show, 
is sufficient to induce an equilibrating process. We have per- 
formed simulations for the singlet decay times after a nuclear 
preparation time r cq (on the order of tens of /is) in which we 
introduce one or two electrons into the double dot with ap- 
propriate bias and magnetic field. Interestingly, in spite of the 
equilibration of the two nuclear configurations, the singlet de- 
cay times are only modestly extended. This is because the 
distribution over If of the collective nuclear states is broad- 
ened during these naturally occurring electron-mediated equi- 



libration processes. The time scale for this nuclear dynamics 
is in the order of hNEz /jf ~ 20 ns, comparable to that for 
the electron spin dynamics, so that quasi-static approximation 
for the nuclear quantities becomes invalid in the regime of 
(Sh z ) « 0. The resulting fluctuations suppress the singlet de- 
cay time enhancement even though Sh z becomes small. In 
contrast, such dynamics does not exist in the fully polarized 
collective states, where the pumping of singlet states pushes 
If to their minimum values. The polarization cycles, while 
unable to produce high degrees of polarization, significantly 
regulate the nuclear spin states and reduce the nuclear field 
fluctuations, thus extend the decay time of the singlet state. 

In summary, we have studied a nuclear polarization scheme 
in gated double dots utilizing the S—T+ degeneracy point, and 
examined its impact on the two-electron spin singlet decay 
time, obtaining two-orders-of-magnitude enhancement. We 
have shown that high degree of nuclear polarization is not es- 
sential to suppress the nuclear relaxation channel for a two- 
electron spin state. Instead, enhancement of the singlet re- 
laxation time is obtained by electron-mediated equilibration 
process within the two nuclear baths that suppresses the Over- 
hauser field fluctuations within each reservoir. We have ex- 
plored other strategies to facilitate this equilibrating process, 
and have found that while equal nuclear polarizations between 
the two dots may be obtained, they are not accompanied by 
narrowing of the nuclear state distribution and thus do not re- 
sult in a dramatic enhancement in the singlet decay time. 

We thank C. M. Marcus and J. R. Petta for very useful dis- 
cussions. This work is supported by NSA/LPS and ARO. 



* Electronic address: guy.ramon@gmail.com 

1 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 
(2004). 

2 A. C. Johnson et al, Nature 435, 925 (2005); R. Hanson et al, 
Phys. Rev. Lett. 94, 196802 (2005). 

3 F. H. L. Koppens et al, Science 309, 1346 (2005). 

4 J. R. Petta et al, Science 309, 2180 (2005). 

5 I. A. Merkulov, Al. L. Efros, and M. Rosen Phys. Rev. B 65, 
205309 (2002). 

6 J. Schliemann, A. V. Khaetskii, and D. Loss, Phys. Rev. B 66, 
245303 (2002). 

7 N. Shenvi, R. de Sousa, and K. B. Whaley, Phys. Rev. B 71, 
224411 (2005). 

8 W. A. Coish and Daniel Loss, Phys. Rev. B 72, 125337 (2005). 

9 C. Deng and X. Hu, Phys. Rev. B 73, 241303(R) (2006). 

10 W. Yao, R. B. Liu, and L. J. Sham, Phys. Rev. B 74, 195301 
(2006). 

11 D. Klauser, W. A. Coish, and Daniel Loss, Phys. Rev. B 73, 
205302 (2006); D. Stepanenko, G. Burkard, G. Giedke, and 
A. Imamoglu, Phys. Rev. Lett. 96, 136401 (2006). 

12 G. Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 59, 2070 
(1999). 

13 J. M. Taylor, C. M. Marcus, and M. D. Lukin, Phys. Rev. Lett. 90, 
206803 (2003). 

14 A. S. Bracker etal, Phus. Rev. Lett. 94, 047402 (2005). 

15 K. Ono and S. Tarucha, Phys. Rev. Lett. 92, 256803 (2004). 



J. R. Petta and C. M. Marcus (unpublished). 

17 R. de Sousa and S. Das Sarma, Phys. Rev. B 68, 115322 (2003). 

18 Zamboni Machine, named after its inventor Frank J. Zamboni, has 
become a generic name for resurfacing ice rink machines. The 
analogy to the supression of quantum fluctuations in the nuclear 
reservoirs was suggested by C. M. Marcus. 

19 X. Hu and S. Das Sarma, Phys. Rev. A 61, 062301 (2000). 

20 Inter-ring dynamics would result in the center rings being po- 
larized first. The outer rings, for which the nuclear dynamics is 
slower, will tend to depolarize the nuclei in the center. Neverthe- 
less, once the polarization is complete, all of the rings' collective 
nuclear states will be polarized, resulting in a qualitatively similar 
distribution to the one that is considered in the main text. Further 
study is needed to quantify the effects of inter-ring dynamics. 

21 C. Deng andX. Hu, Phys. Rev. B 72, 165333 (2005). 

22 J. M. Taylor, J. R. Petta, A. C. Johnson, A. Yacoby, C. M. Marcus, 
and M. D. Lukin. lcond^m at/0602470 (unpublished). 

23 This nuclear state distribution is in fact wider than the distribu- 
tion obtained when considering non-uniform HF couplings, the 
latter given by a sum of Gaussian distributions corresponding to 
the number of nuclei in each ring. 

24 E. A. Laird etal, Phys. Rev. Lett. 97 056801 (2006). 

25 Slicing the dot to 100 rings gives p = 4.3% for N = 10 5 . The HF 
coupling differences within each ring are then smaller than 0.1%, 
roughly preserving the total spin numbers in each ring. 



